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We study the fluctuations of the total current for the partially asymmetric exclusion process in 
the scaling of a weak asymmetry (asymmetry of order the inverse of the size of the system) using 
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the cumulants of the current both the leading and next-to-leading contribution in the size of the 
system. 
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I. INTRODUCTION 

The one dimensional asymmetric simple exclusion process (ASEP) is one of the most simple examples of a classical 
interacting particles system exhibiting a non equilibrium steady state. This stochastic system has been studied much 
in the past, both in the mathematical 0, 0, El and physical literature 0, [H, 0, 0, H, H, EH EH- It consists of particles 
hopping locally on a one dimensional lattice, with an asymmetry between the forward and backward hopping rates. 
The exclusion constraint prevents the particles from moving to a site already occupied by another particle. The 
asymmetry between the hopping rates models the action of an external driving field in the bulk of the system, which 
maintains a permanent macroscopic current in the system. This current breaks the detailed balance and keeps the 
system out of equilibrium even in the stationary state. The special case for which the particles hop forward and 
backward with equal rates is called the symmetric simple exclusion process (SSEP). It corresponds to a situation 
for which the detailed balance holds in the bulk which means (in the absence of boundary conditions breaking the 
forward-backward symmetry) that the system reaches equilibrium in the long time limit. In this case, the system 
belongs to the universality class of the Edwards-Wilkinson (EW) equation 0, E2] ■ On the contrary, if the two hopping 
rates are different, detailed balance is broken and the system reaches in the long time limit a non equilibrium steady 
state characterized by the presence of a current of particles flowing through the system. In that case, the system 
belongs to the universality class of the Kardar-Parisi-Zhang (KPZ) equation [(j [l3j. 

Because of its simplicity, the ASEP is an interesting tool to investigate the general properties of systems out of 
equilibrium. Moreover, the ASEP is related through various mappings to many other models, in particular: the zero 
range process [l4[ , directed polymers in a random medium [l5[ , interface growth models [H, EE EH , the six vertex 
model [TH, E3 , XXZ spin chains (nj EH EH ■ It is also used as a starting point to model physical phenomena such as 
cellular molecular motors [l9l ]. hopping conductivity (2p| . traffic flow [2l| . usually by enriching the dynamics of the 
ASEP by new rules that makes it closer to the studied phenomenon. 

The ASEP, along with a very small number of other statistical mechanics models is known to be "exactly solvable" 
in the sense that several quantities can be calculated exactly, which is a rather uncommon property. The totally 
asymmetric case, for which all the particles hop in only one direction (TASEP), is usually the easiest to solve. On 
the contrary, the model with partial asymmetry often exhibits a more intricate mathematical structure and is thus 
more difficult to solve. A few different app roaches have been used in the past to obtain exact results for the ASEP: 
the matrix product representation p2l . l2.'i| allows to calculate explicitly the probabilities of all the configurations in 
the stationary state for both open systems connected to reservoirs of particles and systems on a ring with periodic 
boundary conditions; techniques from random matrix theory [24], HH, [26[ have been used for the study of infinite 
systems defined on Z; Bethe Ansatz has provided many exact results, principally on a ring 0, but also recently 
for open systems [28j . 

The fact that Bethe Ansatz can be used to study the ASEP is strongly related to the "quantum integrability" 
of the model. Indeed, the Markov matrix governing the time evolution of the probabilities for the configurations of 
the ASEP is similar to the hamiltonian of the XXZ spin chain and closely related to the transfer matrix of the six 
vertex model. While its formulation for the ASEP is well understood, the use of the Bethe Ansatz is usually quite 
technical. The difficulty lies in the determination of the so called "Bethe roots" in terms of which all the quantities 
we want to calculate are expressed. These Bethe roots are solutions of a set of highly coupled polynomial equations 
called the Bethe equations of the system. Their solutions are usually not known for general values of the parameters 
of the model studied. In the case of the ASEP, Bethe Ansatz has been used successfully to calculate the gap of the 
system [l6|, [2|| [29|, [3(J HJ , related to the dynamical exponent which governs the speed at which the system reaches its 
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stationary state. It has also allowed to calculate some properties of the fluctuations of the current in the stationary 

state 0, m, hi m, m, Hi, m . 

In the present work, we study the fluctuations of the steady state current for the ASEP on a ring in the scaling 
of a weak asymmetry between the hopping rates (asymmetry scaling as the inverse of the size of the system). The 
main result of the article is the Bethc Ansatz derivation of all the cumulants of the current in this scaling limit. We 
obtain an explicit expression ([9]) for the leading and next-to-leading contributions in the size of the system. We use a 
rewriting of the Bethe equations for the ASEP in terms of a polynomial functional equation. Our approach is based 
on the functional Bethc Ansatz and does not rely on the behavior of the Bethc roots in the large system size limit. 
We check our results numerically by solving the functional Bethe equation for systems up to size 100. 

In section [IIJ we discuss our formula @ for the cumulants of the current in the scaling of a weak asymmetry. 
In section IIII1 we write the Bethe equations for the ASEP as a polynomial functional equation and recall how this 
equation can be solved perturbatively to calculate the first cumulants of the current. In section HVl we define a new 
version of the functional equation that remains regular in the symmetric limit. Then, in sectionfVl we take the weakly 
asymmetric limit of this equation, and we solve it in section IVII A few technical calculations are relegated to the 
appendix. 



II. CUMULANTS OF THE CURRENT IN THE WEAKLY ASYMMETRIC EXCLUSION PROCESS 



We consider the partially asymmetric simple exclusion process (PASEP) on a ring. It is a stochastic process 
involving classical hard core particles hopping on a one dimensional lattice with periodic boundary conditions. Each 
one of the L sites can be occupied by at most one of the n particles. The system evolves with the following local 
dynamics: in an infinitesimal time interval dt, each particle hops forward with probability pdt and backward with 
probability pxdt if the destination site is empty (exclusion rule) . 



A. Fluctuations of the total current 



We define the total integrated current Y t between time and time t as the total distance covered by all the particles 
in this duration. This is a random variable which depends on the trajectories of the particles starting in some 
configuration C at time and evolving by the markovian dynamics up to time t. We are interested in the fluctuations 
of Y t in the long time limit, when the (finite size) system reaches its unique stationary state which is independent of 
the initial configuration C. We emphasize that we consider here the long time limit for finite systems. We will take 
the large system size limit only in the end. This is a completely different regime from taking the infinite system size 
limit L — ► oo first and studying then the stationary state t — > oo |(| ■ We want to calculate the first cumulants of the 
current with respect to the stationary state probability distribution of Y t , that is its mean value J{x), the diffusion 
constant D(x) and the higher cumulants: 

J(x) = lim VQ- (1) 

t^oo t 

D(x) = lim (2) 

t— too t 

e 3{x) = lim (m^mm±mi . (3) 

t — too t 

The characteristic function, defined as the mean value of e lYt behaves in the long time limit as [37| 

(e^-e-^'*)', ( 4) 

where 7 is the "fugacity" corresponding to the variable Y t . Taking the derivatives of the previous expression with 
respect to 7, we see that -£(7, x) is the exponential generating function of the cumulants of the current in the stationary 
state, i.e. 

t-,/ \ T , x D(x) o E 3 (x) o ,_. 
E(rf,x) = J(x) 7 + ^p7 2 + ^p7 3 + - ( 5 ) 

The generating function of the cumulants £7 (7, x) can be related to the large deviation function of the current G(j, x). 
The function G(j, x) is defined from the asymptotic behavior of the probability distribution of Y t as 

P t (j) = P (— = j\ ~ e - GW ' x )< when t-> 00. (6) 
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The large deviation function of the current G(j, x) is equal to for j = J(x) and is strictly positive otherwise, leading 
to an exponentially vanishing probability in the long time limit except when j is equal to the mean value of the 
current. Writing 

(e"> Yt ) = f dje^Ptd) ~ f dje^- ^^ 1 ~ e *™*i(ry-G(j>)) > (7) 



we observe that E(pf,x) is the Lcgcndrc transform of the large deviations function G(j,x), that is 

E(y > x)=wax(jy-G(j,x)). (8) 
3 



B. Cumulants of the current in the scaling of a weak asymmetry 

The principal result of this article is the calculation of all the cumulants of the stationary state current in the 
weakly asymmetric scaling limit 1 — x ~ 1/L, or equivalently, the Taylor expansion in the vicinity of 7 = of the 
generating function £7(7, x). Using Bethe Ansatz, we find 

1 -E{ 1 x, V )^-e( 1 =^,x = 1-j) (9) 
p p \ L LI 

= 1 + J2 I 2 ^ ~ P ^ M " J V J ' 

In this expression, p is the particle density p = nj L and the function tp is given by 

fc=l v ; 

where the Bj are the Bernoulli numbers. The expression ([9]) gives the leading and next-to- leading order in L of all 
the cumulants of the current by taking the derivative with respect to p. Only the subleading term of © contributes 
to the fc-th cumulant Ef. for k > 3. In this case, E/. is a polynomial of degree k in the rescaled asymmetry v. The 
coefficients of this polynomial are expressed in terms of the Bernoulli numbers multiplied by factorials. The first 
cumulants Ek arc plotted with respect to the rescaled asymmetry v in fig. [T] We note that they show oscillations 
in the parameter v. This oscillation phenomenon of the cumulants of the current has been observed recently in the 
context of electron transport through a quantum dot in (ijj. The special case v = which corresponds to the SSEP 
has already been calculated by Bethe Ansatz in [36|. For arbitrary u, equation ([9]) leads to an expression for the 
large deviation function G up to the order 2 in L, which matches the result obtained in [3(| using the macroscopic 
fluctuation theory developed in [3^, H(| ■ 

We now justify the weakly asymmetric scaling chosen in equation The crossover between the Edwards- 

Wilkinson and the Kardar-Parisi-Zhang behavior lies at a scaling where the asymmetry 1 — x is nonzero but goes to 
zero when the size L of the system goes to infinity. Consider a tagged particle in the system with asymmetry scaling 
as 1 — x ~ 1/L r . During a time interval At, this particle makes a number of rotations R ~ (1 — x)At/L through 
the system. A typical time interval At to consider is the time necessary for the system to reach its stationary state: 
At ~ L z , z being the dynamic exponent of the system. Then, a criterion for separation between weak and strong 
asymmetry is when R ~ 1: when R 3> 1, the system is asymmetric, while R -C 1 corresponds to a symmetric system. 
This leads to 1 — x ~ 1/L Z_1 . The value of z depends on whether the system is in the Edwards-Wilkinson (EW, 
z = 2) or in the Kardar-Parisi-Zhang (KPZ, z = 3/2) universality class. This leads to two natural scalings for the 
asymmetry, 1 — x ~ 1/L and 1 — x ~ 1/ yL. It turns out that both of these scalings are meaningful for the ASEP. 
The scaling 1 — x ~ 1/L corresponds to the weakly asymmetric exclusion process, whereas the scaling 1 — x ~ 1/y/L 
corresponds to the crossover between the EW and the KPZ regimes [I!)]. In particular, the weakly asymmetric scaling 
1 — x ~ 1/L belongs to the EW class like the symmetric exclusion process. 

From the formula ©, we observe that E(p,v) for v ^ is a rather minimal deformation of the case v = 0. This 
modification is in fact needed to ensure that E(p, v) stays invariant under the Gallavotti-Cohen symmetry given by 

Em 



E{j,x) =£(logx-7,x) , 



(11) 
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FIG. 1: Graphs of the first cumulants of the current Ek in the weakly asymmetric scaling limit, obta ined by t aking the successive 
derivatives in p of the generating function ([9]) at p — 0. The rescaled cumulants j^-(ne/ [kL~J p(l — p)) k Eh{v) are plotted 
with respect to the rescaled asymmetry v at particle density p = 1/2 for k from 3 to 11. 



which in the weakly asymmetric scaling leads to 

% I /) = lL-f-^ + ofi),i/j (12) 

for the leading and next-to-leading order. The study of the exact values for the diffusion constant [13, EH and the 
third cumulant [H[ shows that these two cumulants are still given by equation @ as long as 1 — x <C 1/ \f~L. This 
suggests that this minimal deformation is valid in all the Edwards- Wilkinson universality class. 



C. Phase transition 



In [44( , using the macroscopic fluctuation theory for driven diffusive systems developed in [39|, |40( , the existence of 
a nontrivial dynamical phase transition was found in the weakly asymmetric exclusion process, with in particular a 
phase of weaker asymmetry for which the fluctuations of the current are gaussian (at dominant order in 1/L), and a 
phase of stronger asymmetry in which the fluctuations of the current become non gaussian. Let v c be the value of the 
rescaled asymmetry corresponding to the separation between the gaussian and non gaussian phases. For v < v c , the 
large deviation function of the current G(j, v) (respectively the generating function of the cumulants of the current 
E{p, v)) is quadratic in j (resp. /i) at the leading order in the size of the system. Performing the Legendre transform 
of the leading order of the expression ([9]) for E(fi,i/), we find that in the gaussian phase v < v c , the large deviation 
function G(j, v) is given by the quadratic function of j: 



(13) 
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with 

-=p{\-p)v (14) 
P 

-=2p(l-p)L (15) 
P 

at the leading order in L. On the contrary, in the non gaussian phase v > u c , neither G(j, v) nor E(p, v) are expected 
to be quadratic, even at the leading order in L. We emphasize that it does not contradict the fact that the Taylor 
expansion of E(p, v) given in equation @ is quadratic in p at the leading order. It merely means that the large L 
limit of the asymptotic formula for E(p, v) breaks down and does not represent the full function E(p, v) anymore. 
We will come back to this issue at the end of this subsection. 

According to |44] , the density profile adopted by the system is dependent on the value of the current flowing through 
the system. In the gaussian phase, the density profile remains flat for all values of the current. In the non gaussian 
phase however, the density profile depends on the value of the current: there is a critical value j c [y) for the current 
such that if |j| > j c (y), the density profile remains flat, while if \j\ < j c {v), the profile becomes time dependent. The 
signature of this transition between the flat profile and the time dependent profile is visible through the appearance 
of non analyticitics in the large deviation function of the current or, cquivalently, in its Legendre transform, the 
(rescalcd) generating function of the cumulants E(p, v). A non analyticity in the large deviation function G(j, v) at 
j = ±j c (i>) corresponds to a non analyticity in E(p,v) at p Cj i(u) and jJb c ,2{ 1 ') related through the Gallavotti-Cohen 
symmetry as n c ,i(y) + Hc,%(y) = £log(l — vjL) ~ —v. 

We now look at the non analyticitics of the expression ([9]) for E{p, v). From the asymptotic behavior of the Bernoulli 
numbers 



B 2k ~ (-If-HVirk i^—j , (16) 

we observe that ip(z) has a singularity at z = —ir 2 . It corresponds for the function E(p,u) to singularities at the 
points p(l — p)(fi 2 + pv) = —it 2 . This equation has real solutions for fi if p(l — p)v 2 > Air 2 . Thus, non analyticitics 
appears in E(p, v) as soon as v > v c for 

(17) 



and in this case, the non analyticitics of E(p, v) are at the points 



-v±Jv* - 



,2 _ 4tt 2 



»cM») = v 2 P(w) • (18) 

These expressions for p, Ct i(v) and p c ,2(y) should hold only in the vicinity of v = v c since we used the expression @ 
for E(p, v) which is valid for p far from only in the gaussian phase. For the Legendre transform of the gaussian 
leading order of ©, we define the function p{j) such that 

j = L^-E^{j),v) and E(n(j), v) + G(j,u) = j ^ (19) 

By the inverse function of p{j), the values p Ct \(v) and p c ,2{y) correspond for the large deviation function G(j,u) to 
±3c(v) wit; h 



^-*-'V"-^ (20) 

near v = v c . By Legendre transform, the function p(j) sends the region \ j\ < i c {v) (where the density profile is time 
dependent) of the plane (j, v) to the region p Cy i(v) < p < /i c ,2(i / ), or cquivalently p{\ — p){p 2 + \iv) < — ir 2 , of the 
plane (p, v). It gives at the leading order of ([9]) 

L EM < _ n2 
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FIG. 2: Rescaled generating function of the cumulants of the current LE(/j,, v) at half filling, plotted with respect to w = 
fj,/(Llogx) between w = —1.25 and its symmetric through the Gallavotti-Cohen symmetry w = 0.25 for v = 2n, v — An, 
v — 6n, v = 8n and v — Wn. The solid line represents the leading order of the result (|9]) of the Bethe Ansatz calculation. 
The gray dots correspond to the numerical resolution of the functional Bethe equation for n = 25, L = 50. The black dots 
correspond to the numerical resolution of the functional Bethe equation for n — 50, L — 100. 

which agrees with equation (25) of [IH (where p is taken to be equal to 1/2). 

Our resolution of the functional Bethe equation of the ASEP (|2"2"j) in the weakly asymmetric scaling limit, leading 
to the expression ([9]) for the cumulants of the current, relies on a perturbative expansion near fi = of the functional 
Bethe equation. This approach does not always give an information on the value of i?(/x, v) for /i far from 0. 

From the discussion in the beginning of this subsection of the gaussian/non-gaussian phase transition, it is expected 
that the function E{[i,i>) will be equal to its Taylor expansion for v < v c . In particular, the function should be 
quadratic in \i in the large L limit. On the contrary, for v > v Cl E^^v) is expected to be different from its Taylor 
expansion ([9]), even at the leading order in L. 

In order to check whether the function i?(/i, v) was equal to its Taylor expansion in /i = given by ([9]), we studied 
numerically the Bethe equations of the ASEP for systems up to size L = 100 (see appendix |A")) . In fig. [21 we show 
the results we obtained at half filling (y c = 47r) for different values of the asymmetry v and of the size of the system 
L. These results are in excellent agreement with the emergence of non gaussianity for v > An: for < v < An, the 
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numerical evaluation of E(j, x) fits well with the quadratic expression given by the leading order of equation ([9]). For 
v > 47r, however, there is a region between 7 = and its symmetric by the Gallavotti-Cohen symmetry 7 = logs 
where E(p/,x) differs significantly from the quadratic leading order of (|9|). Outside of this region, the numerical 
evaluation of E{-y,x) still agrees with ©. 

We emphasize that the fact E(/i, v) is defined for a finite system in ([5]) and ([9]) as a generating function in that is, 
as a Taylor series for fi at /i = 0, does not contradict the fact that E(/j,, v) can be different from its Taylor expansion at 
/i = in the large L limit. An example of such a behavior is exhibited by the function Ylj=a(~ (1 + I^n 2 ))! / j\ = 

e -1 /^ +0 (1/L). This function of fi and L is, for finite L, a rational fraction in /i which is entirely defined for /i £ C 
by its Taylor expansion in /it = through an analytic continuation, but develops an essential singularity in fj, = 
when L — > 00. 



III. REMINDER OF THE FUNCTIONAL FORMULATION OF THE BETHE EQUATIONS 

In this section, we recall the functional formulation of the Bethc Ansatz for the ASEP. We show how the generating 
function of the cumulants ([5]) can be expressed in terms of a solution to a functional polynomial equation and how 
this functional equation can be solved perturbatively to obtain the first cumulants. 

It can be shown [12, [13] that the generating function -E(7, x) of the cumulants of Yt over the variable 7 is equal to 
the eigenvalue with largest real part of a deformation M^) of the Markov matrix M of the system. Because of the 
underlying integrability of the ASEP, the diagonalization of M(j) can be performed using the Bethe Ansatz. The 
Bethe equations of the system can be rewritten (37| in the functional equation 

Q{t)R(t) = e i7 (l - t) L Q{xt) + (1 - xt) L x n Q{t/x) . (22) 

The polynomial Q is of degree n and the polynomial R is of degree L. We choose the normalization of Q such that 
the coefficient of highest degree of Q(t) is equal to one: 

n 

Q(*) = ![(*-%•) ■ ( 23 ) 

0=1 

If we set t = yi in the functional Bethe equation (|2"2"|) . we obtain 

L 



±ZVi_\ lf^_^. (24 ) 



which is the usual form of the Bethe equations in terms of the Bethe roots y,. Equations (|2"2"1) and (pM|) arc completely 
equivalent forms of the Bethe equations. In particular, they both have a large discrete set of solutions corresponding 
to the different eigenstates of the deformed matrix M (7) . We are only interested in the solution corresponding to the 
largest eigenvalue of M(p/), which is characterized by 

Q(t)=t" + 0( 7 ) . (25) 

Equivalcntly, the Bethe roots yi all tend to when 7 — > for this solution of the Bethe equations. In the following, we 
will also use a relation coming from the fact that the stationary state is a zero momentum state (37j . This condition 
implies 

e"TQ(l) = x n Q(l/x) . (26) 
The generating function E(-/,x) of the cumulants of the current is then given by [37l ] 

E(l,x) n JQ'(1) lQ'(l/x)\_. 1 1 



(,-,) ^.iM-p.^ ' ' . (27) 

V CHI) x u(l/x) J ' V 1 - Vi 1 — xviJ 



p V Q(l) x Q(l/x) J t^V 1 -^ 1 — aryi 

In the rest of this section, we explain how the functional Bethe equation can be solved perturbatively near 7 = 0. 
Introducing the function 



(28) 
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the functional Bcthc equation (|22[) becomes 



x n e (L-2n)j 



e~^R{t) = (1 - xtfA{t) + (1 - t) L . (29) 

Note that we have added an extra factor e™ 7 in the definition compared to [Hj]. In terms of A(t), equations (|25|) . |26|) 
and (|27p become respectively 

4(t) = l + 0( 7 ) (30) 

A(l) = 1 (31) 

M = -(l-x)i'(l), (32) 
P 

In the following, we will also need a few other properties of A(t). Because Q is a polynomial of degree n, we infer 
from the definition ([28]) of that 

lim A{t) = e-™ 7 . (33) 

t — >oo 

This equation fixes the term of highest degree of the polynomial R, using the relation (f29|) between R(t) and 

e' ni R{t) - (a^e - " 7 + x n e (L "" )7 )(-l) L t i is a polynomial in t of degree L — 1, (34) 

Since the value of Alt) is known when 7 is equal to zero (|30[) . it is natural to attempt solving the functional Bcthe 
equation (|29|) perturbatively near 7 = 0. Moreover, a perturbative solution of (|29|) in powers of 7 gives access to the 
first cumulants of the current. Using (|30j) . we write the expansion of A(t) in powers of 7 as 

00 

A(t) = l + ^^.(i)7 fe - (35) 
fc=i 

From (|25|) and the definition l]28[) of we deduce that the Afc(i) are polynomials in 1/i of degree kn. This 

observation will be crucial for the following. We will now sec that equation (|29p can be reformulated as a recurrence 
formula which can be used to calculate explicitly the first Ak(t) and, as a consequence, the first cumulants. In 
particular, the three first cumulants were calculated in [38[ for finite size systems by this method that we now recall 
(see [HI for more details). Reminding that R(t) is a polynomial in t, that is R(t) has only nonnegative powers in t, 
we can eliminate it from equation (|29p by doing the expansion for t — > 0. We have 

Alt) 1 T n„(L-2n)7 
T^TT + 7 U \, ^ = O (t°) . (36) 

{l-t) L {l~xt) L A{xt) v ' v ' 

This equation must be understood in the following way: first, we expand the l.h.s. around 7 = in terms of the Ak(t). 
Then, at each order in 7, the l.h.s. has a finite limit when t — ► 0, that is all the negative powers in t from the Ak{t) 
cancel out. It is crucial to respect the order between the two expansions in 7 and t. Expanding first around 7 = 
makes the poles of A(t) (that is, the Bethe roots yA disappear from the problem, leaving us only with the algebraic 
properties of the polynomials Ak{t). On the contrary, if we did the expansion in t = first, equation (|36p would not 
contain any information since A(t) is regular when t — > 0. Introducing the operator which acts on an arbitrary 
function u as 

(A x u) (t) = u(t) - x ll u{xt) , (37) 
we rewrite the previous equation in the slightly more complicated form 

^(i#y=vM A(it)+ ^) +o( ' r <3s) 

At order k in 7, the r.h.s. depends only on the Aj(t) for j < k. We emphasize that Ak(t) cancels out. This observation 
is the key that allows us to solve A(t) order by order in 7. Noting that A x acts separately on each power of t and 
that A x O (t°) = O (t°), we can invert A„ in (1581. We have 



(l-t) L x \{l-xt) L \ v ' A(xt) J J t 
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The additional term b/t n comes from the fact that the operator gives when applied on \/t n . Recalling (f33|) and 
the fact that the A^{t) have only negative (or zero) powers in t, we finally obtain 



A(t) 



(i - xty 



A(xt) 



; (L-2n) T - 

"Apr. 







- b 


Ul-t) L - 


(t- 1 ) 










(t-°°) 



(40) 



ft 1 

We used the notation [u(i)]/ t _oo) for the negative powers in t of a function u(t) in its expansion near t = 0, after the 
pcrturbative expansion near 7 = as before. The constant b can be set, order by order in 7, by the value pip of 
A(l). We note that in equation (|40|) , L, which was initially the size of the system and the degree of the polynomial 
i?, no longer needs being an integer. It can assume any complex value such that there is a constant b which solves 
A(l) = 1, that is for which the coefficient of b in equation (|40p is nonzero when t = 1. Thus, L must be such that 



n-l 

E 

3=0 



(-1^ = (-1)™- 



L - 1 
n-l 



^0 



(41) 



Wc find that pip can not be ensured if L is an integer between 1 and n — 1, which never happens for the ASEP 
because of the exclusion rule. 

We note that the recurrence (|4"0"|) is singular when x — ► 1 since A x goes to in this limit. More precisely, starting 
from pop , the recurrence equation (|40p gives for an expression which is regular when x — ► 1, as A" 1 is applied 

on an a constant independent of t at this order. Using again (|40p to calculate ^(i), we see that the A" 1 contributes 
a pole of order 1 in x = 1. Iterating (|4"0"|) . we thus see that the Afc(i) have a pole of order k — 1 in x = 1. As we 
arc interested here in a scaling limit for which x goes to 1 as the size of the system goes to infinity, this form of 
the pcrturbative solution will not be usable. We have to transform it to make it suitable for the weakly asymmetric 
scaling limit. 



IV. REGULARIZATION OF THE FUNCTIONAL EQUATION 

In this section, we regularize the function A(t) in the limit x — > 1 and show that the regularized function can be 
calculated perturbatively, in a similar way to what we did for A(t) in a previous section. 

Equations pip and p"!?|) indicate that A(l) is regular in x = 1 and that A'(l) has only a pole of order 1 in x = 1. 
But we also know that at order k in 7 the function A(t) has a pole of order k — 1 in x = 1. This suggests that an 
expansion in t = 1 should allow us to regularize the functional equation (|29p in the limit x — * 1. We define 

A(y) = A(l - (1 - a:)i/) . (42) 

In appendix [B] we argue that A(y) is regular when x — * 1 and work out explicitly how the cancellation of divergent 
terms in x = 1 works in the case of a system with only one particle. In the rest of this section, we will write a 
recurrence equation for A(y) similar to (|4"0"]) for A(t) but regular in the limit x — * 1. This will allow us to study the 
weakly asymmetric scaling limit. A difference will be that we will now have to do expansion both in 7 and 1 — x. 

A. Recurrence equation for A(y) 

In terms of A(y), the functional Bethe equation (|2"9"|) rewrites 

e-"TR(l - (1 - x)y) ,„ T , . , x n e^- 2n ^ 

(1-xT A(l + xy) 

while equations p0|) . (J3TJ) and (|32|) become respectively 

%) = 1 + 0(7) (44) 
1(0) = 1 (45) 

fefl=l'(0). (46) 
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We will now solve the functional equation (f43|) perturbatively near 7 = and x = 1. We first write the expansion of 
A(y) near 7 = and x = 1 as 



(47) 



fc=l z=o 



The range for the summation over k comes from (|44j) while the range for the summation over I is a consequence of 
the regularity in x = 1 of Since the Ak(t) are polynomials (of degree kn) in 1/i, we see that A(y) has only 

positive integer powers in y after the expansion near 7 = and x = 1. A crucial point in the following will be that 
the Ak,i{y) are in fact polynomials in y of degree k + 1 — 1, as can be seen writing Afc(l — (1 — at order Z in 1 — 2; 
in the following way: 



iijfc,i(l/) = [i4*(l-(l-x)y)] (1 _ x) i 



£ ( -iy ( i-z)V^#(i 

J=0 



(i-x) 1 



k+1-1 



(1) 



(48) 



(l-x)'-J 



We used the fact that the Af~ (t) have a pole of order k — 1 in x = 1 . We will now eliminate the polynomial R from 
the functional equation (|43[) in a similar way to what we did in the previous section for the functional equation (1291) . 
We divide the functional equation (|4"5|) by y L (l + xy) L and make the expansion y — > 00. Taking (f3~4"|) into account, 
we obtain 



1 



„n (i-2n)7 



^(y) 

y L (1 +xy) L i(l + X y) 



(e- 



x n-L e (L-n)i\ 



= O 



1 



(49) 



which must be understood as: each term of the expansion in power series in 7 and 1 — x of the l.h.s. is of order \/y L+1 
when y — * 00. Once again, we see that the polynomial R has disappeared, leaving us with an equation involving 
only A(y). Replacing y by y/x and dividing everything by x n , equation (|49j) can be rewritten in the slightly more 
complicated form 



A(y + 1) A(y) = x L - n A(y/x) - A(y) 
(y + 1) L y L y L 



1 



(y + 1) 1 



A(y+1) + 



3 (L-2n) 7 



A(y + 1) 



— x L ~ n e~ n ~ < — e ( - L_ ™' )7 



O 



,L+1 



Defining the finite difference operator A acting on an arbitrary function u as 

(Au)(y) = u(y + l) - u{y) , 

the functional equation (|50|) finally becomes 



My) 



U(y) | V(y+1) 
y L (y + 1) L 



O 



1 



,L + 1 



with 



and 



U(y) = x L - n A(y/x) - A(y) 



V(y) = A{y) 



; (i-2n) 7 

My) 



(50) 
(51) 

(52) 

(53) 
(54) 



Similarly to what happened in the recurrence equation (|38[) for A(t), at order k in 7 and Z in 1 — x, the r.h.s. of (|52[) 
depends only on the with either « = k and j < Z (U(y)) or i < k and j < Z (y(y)). Thus, equation (|52|) provides 
a solution order by order in 7 and 1 — x of ^4(2/)- 



B. Inversion of the operator A 

We see that contrary to what happened in (f3"8"|) . the operator A acting on A in the l.h.s. of (fB"2"l) docs not vanish 
when x — > 1. This is related to the fact that that A(jj) is not singular in the limit x — > 1. The operator A acts 
formally as 

A = - 1 with Dj, = d/dy . (55) 
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Using the Taylor expansion 



D 



1 ^ 7 



3 -z* 



(56) 



where the Bj are the Bernoulli numbers, we see that we can invert the operator A in f|52|) by multiplying both sides 



by D y /(e D « - 1). We obtain 



D u 



My) 



(u(y)\ 
^ - i V y L J 



-D 



v (v(y)\ 



e °y — 1 \ y L J 



O 



,L + 1 



(57) 



The differential operators D y /(e Dy — 1) and —D y /(e~ D y — 1) must be understood as the formal series (|B"6"|) in D y . 
We now define some notations that will be useful in the following. For any function /, we write the series expansion 
of f(y) when y — > oo as 



oo 1 
f(y) = ^2lf}(r) — 



(58) 



where a is a (possibly negative) integer. We will note [/(y)](-) the singular part of f(y) when y — > oo and [f(y)](+) 
the non-singular part of f(y) when y — * oo, that is 



[/(v)k-)= Etftai 



(59) 



and 



[/(»)] 



OO _^ 



(60) 



r=0 



When the function / depends also on 7 or x, we define [/(y)](_) and [/(?/)](+) such that all the expansions in powers 
of l/y when y goes to infinity must be done after the expansion in powers of 7 and 1 — x. Using (|47p . we can expand 
U(y) and V(y) near 7 = and x = 1. Since the Afc^y) are polynomials in y, we see that at each order in 7 and 1 — x, 
U(y) and V(y) are also polynomials in y. We can write 



(r)I/ 



(61) 



r=0 



r=0 



With these notations, we integrate equation ([5T|) as a formal series in y and only keep the divergent powers when 
y — > 00 (that is, the strictly positive powers in y). Taking (|45|) into account, we obtain 



A(y)-i = jt\y L I dy([u] (r) 



D 



V y r-L 



(62) 



(-) 



We did not add a constant term when doing the integration. This can be understood using the analytic continuation 
for complex L in our formulas: since the Ak,i(y) are polynomials in y, they have only positive integer powers in y. 
But adding a constant term when integrating (|57[) would contribute a y L to A(y). This is not possible as we have 
seen previously that we could take for L any complex value. Using (|56j) . we rewrite (|62(1 as 



00 00 



i=0 r=0 



r- L 
J 



y L / dy{[U] [r) + {-iy[V] [r) )y 



r—L—j 



(63) 



(-) 



After calculating the integral and using 



KTH-^tT 1 )' 



we find 



L + j - r - 1 

3 



r=0 j=0 



L + j- r 



^((-ly^W + Mw) 



(64) 
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We observe that the denominator L + j — r — 1 in this equation makes the previous equation divergent if L is taken 
to be an integer. Thus, equation (|64[) must be understood by an analytic continuation in L. The binomial coefficient 
in (|64|) is thus a polynomial of degree j in L. For j ^ 0, we see that the denominator L + j — r — 1 cancels with 
the binomial coefficient, giving non singular terms in the limit where L becomes an integer. For j = though, the 
denominator L + j — r — 1 does not cancel with the binomial coefficient, ft seems to make it impossible to have a 
finite limit for the r.h.s. of equation (|64|) when L tends to an integer. But, as we know that A(y) is analytic in L 
provided that L > n, it simply means that the numerator contains factors canceling the non-analyticities in L integer 
larger than n. This is indeed shown by using the fact that the Ak,i(y) are polynomials in y of degree k + I — 1. From 
the definitions ([53)) and (|54|) . this implies that U(y) and V(y) are polynomials of degree fc + I — 2 at order k in 7 and 
I in 1 — x. Thus, only the terms with r < fc + I — 2 contribute to ;(y) in equation (|64p . If we choose k and Z such 
that k + I < L, we will only need the terms of the sum over r such that r < L — 2. We observe that for these terms, 
the denominator L + j — r — 1 is always nonzero, even for L integer. Thus, using the notation 



L+l 



O (7, 1 - x) L+1 ^J2° (^) ° H 1 x) L+1 ~ S ) . (65) 



s=0 



we can rewrite (|M|) as 



L-2 



L + j - r - 1 

%) -! = -££ Bj V L+ 7^—^- y r+1 - j ((-iy[U] {r) + [V] {r) ) + O ( 7 , 1 - . (66) 



L + j — r — 

r=0 j=0 ■' 

In this last equation, we can take L to be an integer: there are no divergences anymore. 



V. RECURRENCE RELATIONS IN THE WEAKLY ASYMMETRIC SCALING LIMIT 

In this section, we take the weakly asymmetric limit [x = 1 — v/L, L — > 00) of the recurrence relation (|66|) for 
A(y). We write closed equations verified by the leading and next-to-leading expressions in the size of the system of a 
rescaled version of A(y). 

From the expression (|54[) for V(y), we see that we will have to take 7 of order 1/L in order to obtain a non trivial 
expression at finite density p = n/L. It is useful to consider the function h(y) 

m . Kk „, . = = (67) 

From now on, we will no longer use the variables 7 and x. All the expansions in powers of 7 and 1 — x will be replaced 
with expansions in powers of the rescaled variables fi and v. Recalling the fact that the Akj(y) are polynomials in y 
of degree k + 1 — 1, we obtain from (|4"T|) that h(y) is of order 1/L when the size of the system goes to infinity. We 
write 

%) = ^ + ^WtV) ■ (68) 



h(») = 0(M) (69) 
h(0) = , (70) 



Equations (|44[) and (|45|) become 



while from (|46|) . the rescaled generating function of the cumulants of the current E{p,v) is given in terms of the 
function h by 

^ U) ^^ 1 -E(,= ^ X = l-^=h'(0). (71) 



p L^cop V L L 
The binomial coefficient appearing in (|66|) has the L — > 00 expansion 

L + j — r — 1 
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where l r <L_2 is equal to 1 if r < L — 2 and otherwise. Inserting (|72| in (j66|) . the recurrence becomes 



%> = -j 



r=0 



J=0 



1 /y 2 d 2 



L\2 dy 2 



dy 



(-!)'[£%) + M(r; 



(73) 



(-) 



Using (|56p to perform the summation over j, we find 



r=0 



(r) 



(r) 



- i e Uv - i 



J (-) 



We can finally perform the summation over r and we obtain 



Kv) 



i 



U(Ly) V(Ly) 



e-Vv - 1 e 1 /* - 1 



(74) 



(75) 



1 

I 2 
1 

I 2 



2 \dy 2 e~ x lv - 1 



2 \dy 2 e l /y - 1 



(-) 

U(Ly) + y 

V(Ly)+y 



1 



dy e-Vw - 1 J dy 
d 1 



(yf/(iy)) 



(-) 



e 1 /^ -l) dy 



-r{yV(Ly)) 



(-) 



We recall that, from (|53[) . (j54")l and (j6"T|) , U(Ly) and V(Ly) are well defined functions depending on y only through 
ft(y). Equation (|75|) is thus a closed equation verified by /i(y). It will allow us to obtain both the leading (ho(y)) and 
next-todeading {h\{y)) terms of h(y) in the size of the system. We now have to expand both U(Ly) and V(Ly) to 
order two in 1/L to obtain an equation for ho and an equation for hi. The definition of U (|53|) gives 



U(Ly) = U^ {y) + E^M±I^ 



O 



(76) 



with 



{y) = {e^ 1 ^ + h Q {y)) 
U (11) {y) = {e-^- p > - 1) 

U^\y) = e-^"> (uyh' (y) - ^^(1 + 
where p = n/L is the particle density. Similarly, using the definition of V (|54j) . we find 

V(Ly) = X[Lh(y)] - - e™» = V^(y) + K^M^) + V<"%) 



(77) 



where we have defined 



and 



X[f] = 1 + / 



; (l-2p)/i 



(78) 



(79) 



V (0 \y) = X[h (y)} - e -PK- e -0-P> _ e (i"P)M 

V< u >(y)=x'[fto(i0] 



(80) 



ydo) (2/ ) 



From (|77|) and ([80]) . equation (|75|) gives at the leading order in the following closed equation for ho(y): 

- e i/y f /(o)( y ) + y(o)( y) - 



ho(y) 



A/V - 1 



(81) 
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In the next section, we will solve this equation to find ho(y) and thus the leading order of the rescaled generating 
function of the cumulants of the current. At the next-to-leading order 1/L 2 , equation (|75p leads to 

[(l-a(y))/n(y)](-) = [/3(y)](-), (82) 

with 
and 

2\dy 2 e 1 /y~l) W U \dye l /y-l)dy yy W>> 

We used the fact that h\{y) = since h(Q) = ((70)) . Once we know ho(y) by solving equation (f8"Tj) . equation 

(|82|) becomes a closed equation for h\{y). 



VI. SOLUTION OF THE WEAKLY ASYMMETRIC FUNCTIONAL RELATIONS 

In this section, we explicitly solve the functional relations (f8"Tj) and (f8"2")) to all order in [i and v. We obtain the 
leading and the next-to- leading order in L of the rescaled generating function E(/j,, v). 

A. Leading order 

Inserting the expressions (|77|) . (|75|) and (|5D|) for U and 1/, and using [ho(y)],\ = ho(y), the equation (|8ip for the 
leading order becomes 

P (l-2p)M \ 1 p -pn p -{\-p)v , (l-p)p / 1 \ 

.*.-*-*•(! + M»» + ^) = e el/ ,_7 + ° (£) ■ (85) 

Multiplying both sides of the previous equation by y ^e^» — e ^ =1 + we obtain 

ye ^ e -(i-p)-(l + hQ{y )) + I l + ho{y) = V (e~^e-^-e> + e^^j + O [^o) ■ ( 86 ) 
Let us write the unknown function ho(y) as 

, . . „ (l-2p) M + (l-p)^ r(y)-l 

h (y) = -l + e s e 2 « . (87) 
We can state ([56]) as an equation for r(y). We obtain 

>+ (1 - 



, r(y) 

y cosh 

2y 



= j/ cosh 

(-) 



Because of (|70j). we also want ho(y) to have only strictly positive powers in y order by order in [i and v. This can be 
expressed in the form 

[ho(y)}(+) = . (89) 

In terms of r(y), it becomes 

_ (l-2p)M-(l-p)„ 

= e 2 . (90) 

(+) 



e 2 » 



We prove in appendix [Cl that 

r{y) = y/1 - 2(1 - 2p) W - 2(1 - p)yv + y 2 ([i + (1 - p)^) 2 (91) 

is the unique formal power series in \x and v that solves (|88[) and (|90p . Thus, equation (|87|) for ho[y) with r(y) defined 
by (|91[) is the unique solution of (|8"Tj) . 
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B. Next-to-leading order 



With the expression (f87|) for ho(y), we can simplify the expressions ([83]) and (|84[) for a{y) and f3(y). We begin with 
a(y). Inserting in (|83|1 the expressions i|77]) and dHHI) of E/ (11) (y) and 

we obtain 



a(y) 



-(i-P> _ i X'[/i (y)] 



e-Vv - 1 eVi/ - l 
The expressions of the operator X (179)) and of ho(y) l[87|) give 

1 _ e - r {y)/y 



l-a(y) = e-^- p >- 



i - e-yy 



We simplify in appendix [Pi We find 



(l-p)v 2 e 



1p-pii p -0--p)v 



4 sinh 



(4) 



(92) 



(93) 



(94) 



Mi 



Ay 2 r(y) 



sinh 



(*) 



To be able to solve ([82]) . we will need to factorize 1 — a(y) as a product of a function u(y) with only positive (or zero) 
powers in y and a function i>(y) with only negative (or zero) powers in y, after the expansion in powers of \x and v. 
A possible factorization is 



and 



u(y) = e- {1 - p}u r(y)e^ 
sinh 



r(y)sinh (i 



(95) 



(96) 



After the expansion in powers of /i and ^ , r(y) is equal to 1 plus strictly positive powers in y (|9ip . which shows that 
u(y) has indeed only nonnegative powers in y. On the contrary, the fact that the hyperbolic sinus has only odd powers 
gives 



sinh 



OO ^ 

S (2.7 + 1) 



r(y)/2y ^ (2j + 1)! V% 2 



1 (l-2p)/i (l-p> , (/i+(l-p)i/) 



2 \ 3 



2y 



2y 



(97) 



Together with 2ysinh(l/2y) = 1 + (1/y), this proves that w(j/) has only negative (or zero) powers in y. Writing the 
equation ([82]) for the next-to- leading order of h(y) as 



(l-a(y))h 1 (y)=/3(y)+0 



and using the factorization of 1 — a(y), we divide the previous equation by v(y). We obtain 

u(y)h 1 (y) = ^\ + 
v{y) 

Noting that u(y)hi(y) has only strictly positive powers in y, we can write 

"/%) 



Dividing by u(y), we finally find an expression for h\{y) 

hi(y) 



v{y) 



(-) 



u(y) 



v (y) 



(98) 



(99) 



(100) 



(101) 



(-) 
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C. Calculation of the generating function of the cumulants of the current 

We write the rescaled generating function of the cumulants of the current as 

Ejjh v) = EiQ^v) E 2 (p,,v) 
p L L 2 

Since equation (|7ip expresses E(/j,, v) in terms of the derivative of h(y) in y = 0, we need the expansion of h{y) at 
first order near y = 0. We also note that the expression pOip for h\{y) involves the singular part (positive powers in 
y) of the expansion in l/y when y — > oo of a function of y. However, this is not a problem: as h(y) is a polynomial in 
y at each order in fi and v, its expansion when y — > oo has a finite number of terms, which are all positive powers in 
y. Using the value ([57)1 of ho(y), the generating function (fTTj) at the leading order in the size of the system becomes 

Ei(/i,u) = P (l - p)(fi 2 + H . (103) 
For the next-to- leading order, we need the expression (|101|) of hi(y). Using (f95|) , we have 

(1— 2p)u— (1— p)v 

u{y)=e- V +0(y). (104) 

The next-to-leading order of the eigenvalue is then 

(105) 

with the notation [f(y)](i/ y )° for the constant term in the expansion in l/y, after the expansion in powers of /z and 
v as usual. The calculation of this expression is done in appendix [E] Wc find for the next-to-leading order of the 
eigenvalue 

& ( \ p(l~p)M 2jy . B 2 k-2 kn vfc/ 2 . \k nw\ 

E 2 (n,v) = + 2^ fciffc- 1V P ^~ p ' & + ^ v > ' ( 106 ) 

fc=i 

which concludes the proof of equation j9]) . 



(102) 



Ei(p,,v) 



-(l-2p)n + (l-p)y 



yv(y)_ 



VII. CONCLUSION 



Exact results for the cumulants of the stead y st ate current in the exclusion process on a ring have already been 
obtained in the past using Bethe Ansatz: in [361 ] . all the cumulants have been calculated in the thermodynamic 
limit for the symmetric exclusion process, while in [37l. [381] . finite size expressions were obtained for the three first 
cumulants in the system with partial asymmetry. In this paper, we calculated all the cumulants of the current when 
the asymmetry scales as the inverse of the size of the system (weakly asymmetric exclusion process). Wc obtain for 
all the cumulants both the leading and next-to- leading contributions in the size of the system ^ . 

In the scaling of a weak asymmetry, it has been pointed out recently [44j that the system exhibited a non trivial 
phase diagram, with in particular a phase of weaker asymmetry in which the fluctuations of the current are gaussian, 
and a phase of stronger asymmetry for which the fluctuations become non gaussian. On our exact formula @ for 
the cumulants of the current, we observe that the next-to- leading order develops singularities when the rescaled 
asymmetry v is larger than some critical value v c , in perfect agreement with what was predicted in |44| on the basis 
of the macroscopic fluctuation theory [3^, H(| which provides a hydrodynamic description for a large class of driven 
diffusive systems. Moreover, from a numerical solution of the functional Bethe equations for systems up to size 100, 
we confirm that the fluctuations of the current become non gaussian if the asymmetry parameter is larger than v c . 
This fact can unfortunately not be seen on the exact formula ([9]) for the generating function of the cumulants as the 
non gaussianity of the fluctuations of the current at the leading order is not encoded directly in the large system size 
limit of the cumulants of the current but is hidden in the non perturbative behavior of their generating function. It 
would be interesting to calculate by Bethe Ansatz the full form of the generating function of the cumulants, including 
its non perturbative behavior. 

We observe on the generating function of the cumulants ^ that the weakly asymmetric case is given by a small 
deformation of the generating function of the symmetric case, even if the asymmetry parameter becomes larger 
than the critical value v c . This deformation can be understood as a minimal way to preserve the Gallavotti-Cohen 
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symmetry. If we go further away from the symmetric case, it is known that the cumulants of the current have more 
complicated expressions. When the asymmetry scales as the inverse of the square root of the size of the system 
(crossover between the Edwards- Wilkinson and Kardar-Parisi-Zhang universality classes), the three first cumulants 
are indeed given by multiple integrals [H, [43| . It is still an open question to calculate all the cumulants of the current 
in this crossover scaling. 

Our method for solving the Bethe equations of the system is different from the one used in [36] for the symmetric 
case. In that article, the authors used directly the expression of the Bethe equations in terms of the Bethe roots. 
Their method relies on the fact that the Bethe roots accumulate on a curve in the complex plane as the size of the 
system goes to infinity. For general Bethe equations, it is in general difficult to know what this curve is: it usually 
requires a numerical resolution of the Bethe equations, which is not always possible since the Bethe equations are 
highly coupled. The method we use for solving the Bethe equations, in contrast to the one used in [301, does not 
rely on the behavior of the Bethe roots in the large system size limit. Instead, we use the formulation of the Bethe 
equations as a functional polynomial equation known as Baxter's TQ equation. This equation can be solved, in the 
case we are studying, by purely algebraic manipulations. It would be interesting to know if such an approach could be 
used to study the Bethe equations for some other problems. In particular to calculate the fluctuations of the current 
for the case of the open ASEP [28| and the multispecies ASEP [45( for which the Bethe equations are already known. 
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APPENDIX A: NUMERICAL SOLUTION OF THE FUNCTIONAL BETHE EQUATION 

We used Newton's method to solve the functional Bethe equation starting from the known solution at \i = 0. 
However, because the coefficients of the polynomial Q(t) do not vary slowly with respect to 7, we were not able to 
perform our numerical study on the original functional Bethe equation (|22[) . Instead, we used an equivalent equation 
which does not involve Q(t) but only R(t) [46[. This equation can be obtained from (f2"2"| in the following way: first, 
we divide (|2"2")l by Q(t) and obtain 

R(() = ^ (1 -^ + <i-*0V.M£>. (Ai) 

Then, we replace t in the previous equation by t/x. We have 

fl(fW =^(l- f A)^ + (l-*)^§^. (A2) 

Multiplying (|A1[) and (|A2[) . we see that among the four terms coming in the r.h.s., only one is not proportional to 
(1 — t) L . Moreover, all the Q(t) and Q(t/x) cancel in this term. This leaves us with the equation 

R{t)R{t/x) = x n e L ^(l - xt) L (l - t/x) L + O ((1 - t) L ) . (A3) 

This equation provides L constraints on the polynomial R. Adding the additional equation 

R(0) = x n + e i7 , (A4) 

which is a consequence of (|22[) . or 

R(l) = e" 7 (l - x) L , (A5) 

which is a consequence of (|26|) and (|22|) . we have enough equations to constrain completely the polynomial R(t), 
which is of degree L. It turns out that the coefficients of R(t) vary much more slowly than the coefficients of Q(t), 
allowing us to perform a numerical calculation of E(-f, x) for systems with up to 50 particles on 100 sites. We did these 
calculations keeping 400 significant digits for the coefficients of R(t). In the end, we obtained a numerical evaluation 
of E(^/,x) which was symmetric through the Gallavotti-Cohen symmetry, which validated our numerical calculation. 
A plot of the result for E(^,x) is shown in fig. O The results are discussed in section HT1 
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APPENDIX B: REGULARITY OF A{y) IN x = 1 

In this appendix, we explain why A(y), defined in (|42p . is regular in x — 1. In the algebraic formulation of the 
Bethe Ansatz (see [10|). one defines a transfer matrix r(A) which commutes with the Markov matrix for all complex 
value of the spectral parameter A. This transfer matrix can be seen as the generating function over the variable A 
of non local (and non hermitian) generalized quantum hamiltonians (47[. The Markov matrix is given in terms of 
the transfer matrix by the relation M = pe 7 r'(0)r _1 (0). The largest eigenvalue of the transfer matrix e(A) can be 
expressed in terms of the Bethe roots j/j. An expression for e(A) is given e.g. in [Icj . equation (68), for 7 = and in 
terms of the Bethe roots = e 7 (l — yi)/ (1 — xy{) (the authors use the convention p + q = 1 and p and q are exchanged 
in comparison to our notations; the generalization to 7 7^ is straightforward). Defining the variable t in terms of 
the spectral parameter A as 

t= 1 - e " 7 * (Bl) 
we find that e(A) can be expressed in terms of Q{t) and R(t) as 

K ! \l-xtj Q{t) Q(t) {l-xt) L ' v 1 

In the same way that the largest eigenvalue of the deformed Markov matrix is well defined for x = 1, e(A) is not 
singular for x = 1. In particular, its successive derivatives in A = correspond to the largest eigenvalue of a well 
defined generalized hamiltonian and must be regular at x = 1. In terms of A(y), the eigenvalue e(A) rewrites 

e(X)=A(y) + X L ^- -, (B3) 

A(l + xy) 

while the spectral parameter A can be expressed in terms of y as 

1 + xy 

Taking the successive derivatives of e(\(y)) with respect to y, we have 

^(y)) = E^} (") n X ^vZ . (B6) 

dy ^O 1 ) 1 \ J J (l+xy)3+ k d\J w ' v J 

which can be checked by recursion on k. Inserting the expression (|B3[) for e(A(y)) in the previous equation, we observe 
that the term with A L /yl(l + xy) does not contribute for y = if k < L. We obtain for the fc-th derivative of A{y) at 
y = (k < L) 

This shows that A(y) is regular in the vicinity of x — 1 like e(A), at least up to order L — 1 in y. To confirm this, we 
calculated A(t) and A(y) using equation for all systems of size L < 15 and ?? < L/2 up to order 8 in 7. In all 
these cases, we verified that A(y) is regular near x = 1 at all order in y. In the rest of this subsection, we write the 
complete expressions of Q(t). A(t) and A(y) for systems with one particle on a lattice of size L. Again, we observe 
that the expansion of A(t) in powers of 7 is singular in x = 1 while the expansion of A(y) is regular. For n = 1, Q(t) 
is a polynomial of degree 1. It can be written as Q(t) = t + Q(0), the constant Q(0) being set using |26|) . We find 

Q(t) = t-^-, (B7) 

which does not depend on the size of the system: a particle only feels the finiteness of the lattice through its interactions 
with the other particles. The generating function of the cumulants of the current is given (|2T[) by 

E(j,x) = {l-e- 1 )^ ~x) . (B8) 
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Using (28]), we find for A(t) 



X (a z-e^) 



(1 - x)t - x(l - t)(l - e" 7 ) 



^(*) = V : ./n = \i i , Lv 1/ ■ ( B9 ) 



eT t- 



(l-x)t-(l-i)(eT-l) 



The function Ait) has one pole and one zero, which both tend to zero when 7 — > 0. The beginning of its expansion 
near 7 = is 

- 4 « = 1+ (- 1+ 7+ G + 5(^)i + a4?) 72 (B10) 



1 7-2x + x 2 2 1 

6 + 6(l-z) 2 t ~ (1 - x)H 2 + (1 - x)H 3 



7 3 + O ( 7 4 ) 



This expansion in singular when x — ► 1. Taking f = 1 — (1 — x)y in the expression for -A(i), a factor 1 — 1 cancels 
between the numerator an the denominator, and we obtain for A(y) 

1 - (1 - x)y - y(e-r - 1) 1 ' 

Making the expansion in powers of 7 and 1 — x, we find 

A(y) = 1 + ((1 - x)y + (1 - x) V + (1 - xfy 3 + O ((1 - z) 4 )) 7 (B12) 
+ (» + (1 - x) (-§ + 2y 2 ) + (1 - z) 2 + 3y 3 ) + (1 - *) 3 + 4y 4 ) + O ((1 - *) 4 )) 7 2 

+ (V + (1 - X) (| + 3y 3 ) + (1 - X? + 6 y^ + (1 - S )3 ^ + lOy^ + O ((1 _ X )*)\ 7 3 + O ( 7 4) . 

This expression is indeed regular when x — > 1. 



APPENDIX C: PROOF OF THE EXPRESSION ([9TJ) FOR r ( j/) 

In this appendix, we prove that the expression (|9ip for r(j/) is the unique formal series in /1 and v which solves 
equations (f88|) and ([90]) , and such that (|69|) holds. We first check that r(y) given by (j9Tj) solves equation (|88|) by 
direct substitution. We have 



, r (y) 

y cosh — — = y . 

2y y j^{2j)\\2yJ 



^J_(rhj)\ 2j = 



3=0 



( M +(1-P)^) 2 (l-2p) M +(l-p) I / | 1 ^ 
4 2y 4y 2 



(CI) 



The cosh has eliminated the square root of r(y). Taking now only the strictly positive powers in y, we obtain 



y cosh 



r(y) 
2y 



1 /( M +(i-p» 



(-) 



2\ 3 



y cosh 



/i + (l-/3)v 



(C2) 



which is equation (88)) . For equation (|90p . we must as usual do the expansion of r{y) in powers of n and v before the 
expansion in powers of 1/y. Thus, we must write 



r{y) - 1 = (l-2 P )y+(l- P )v 
2y 2 



(O (/i 2 ) + (v 2 ) + 0(^)0 (1/)) £ £ 2/i\*(2/)/iV 



(C3) 



fc=0 z=o 



where the Pk,i{y) are polynomials in y. Taking the exponential of the last equation and expanding again in powers of 
fi and v the term with the double sum over k and I, we obtain 



e 2 » 



(l-2p) M + (l-p)i; 

: e 2 



1 + (O ( M 2 ) + O (v 2 ) +0(n)0 (1/)) £ £ yPkMW 



(C4) 



/c=o ;=o 
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Taking the nonpositivc powers in y only eliminates the term with the double sum, leaving us with equation (|90|) . 

The unicity is obtained from the following argument: let us assume that there exists another solution s(y) of the 
equations ([88]) and (|90|) such that h (y) = O (fi). Because of ([87) . the condition h (y) = O (fi) gives 



s(y) = l + 0(p)+0{v) 



(C5) 



while equation ([88]) gives 



y cosh 



r(y) 



2y 



(-) 



y cosh 



g(y) 
2y 



(C6) 



(-) 



and equation (|90|) gives 



r 




r i 


e 2 » 


(+) 


e 2 « 



(+) 



(C7) 



Expanding the last two equations at power fc in [i and Z in z/, we find by recurrence on k and / that [r(y)](_) = [s(y)h_) 
and [r(y)/y}(+) = [s(y)/y](+) at all order in fi and v. Thus, r(y) and s(y) are equal at each order in /i and v which 
proves unicity. 



APPENDIX D: CALCULATION OF j3{y) 



In this appendix, we simplify the expression ([B"4"[) for f3(y), taking into account the value (|57|) of ho(y). Using the 
identity 



y 2 ( d 2 



i 



2 Vdy 2 e* 1 /* - 1 



f(y)+y 



l 



dy e ±:L /i/ - 1 J dy 



Uyfiv)) - v(d2 



y \ tt..\ , 2 ( d 1 



2 Vrfy 2 e±V» - 1 



f(y)+y 2 



\dy e ±x /v - 1 



f(y) (Di) 



valid for an arbitrary function /, and the expressions ([77)1 and ([80]) of U^°'(y), £/( 10 )(y), V^°\y) and V^ 10 \y), we 
rewrite (1841) as 



P(y) 



-^~p> {vyh'M - + ho(»))) (i-p)^ e - We -(i- rt ^ 



- i 

^M 2 



- 1 



y ( j£ y 

2 \dy 2 e 



2 \dy 2 e 



(I/) 



e -pM e -(l-pV _ 6 (1-P)^ 



^dy eVy -lj dy 

We use the definition ([79]) of the operator X to express the equation ([81]) for /io(y) as 

y*[/k>(tf)] = y (e (1 - p) » + r^-C-*) + y (l - e^ 1 "")^) (1 + /> (y)) + ° (jo ) > 



(D2) 



(D3) 



and its derivative 



d 
dy 



! _ e -(i-p) " e i/» 



(i + fto(y)) 



(D4) 



These last two equations allow us to eliminate all the X operators in the expression (jD2[) of /3(y). We obtain 



fry) = -^Lf>l- e -P^-P> 1- (1-P)j/V (1 _ p) v ±iyL e -(l-rt^l/«^ (y ) +0 f IV 

2 ( e i/a _ i) 2y 2 (e 1 ^ - 1) ^ V " e 1 ^ - 1 ou " / 

1 (D5) 

From the explicit expressions (|8T[> and (|9 If) for ho(y) and r(y), the function h' (y) can be written in terms of 1 + ho(y) 
as 



h , f , l + h (y) 1 - (1 - 2p)y^ - (1 - p)yv 

K{v) = V Ky) 



(D6) 
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We insert this in equation (|D5[) . It leads to 

( v{\-{\-p)yv) (\-yv){\-(\-2p)yp-(\-p)y V ) \ e -(i-^ e i/« / 1 

" I ^ + V^) J eW-i (1 + My)) + ° U 

The first term of f3(y) in the last equation is equal to — y(l — p)v 2 e~ PIJ, e~( x ~ p > v /2 + O However, it is better 

at this point to write all the terms of /?(y) with 1/ sinh(l/2y) in factor. Noting that 



="•-1 2smh(i) 2smh(i) V»». 

and using the expression (|87|l for 1 + ho(y) in terms of r(y), we finally obtain the following expression for (3{y) 
(1 - o)v 2 e~ pp e~ { - 1 ~ p ' )v 

m = (1 P )ve e (Dg) 

4Binh(^) 

- (1 - p)yi/) , (1 - yv){\ - (1 - 2p)y/x - (1 - ^ e (1 - 2p) V <1_p) " e^f 1 / 1_ 



% 4y 2 r( 2/ ) ^ s "i h (^) 



APPENDIX E: CALCULATION OF E 2 (p,u) 

In this appendix, we calculate the next-to-leading order of the generating function of the cumulants in the weakly 
asymmetric scaling, starting from (|105p . Using (|94|) and |96|) . we have 



{l-p)v 2 e M+<1 2 P> r{y) v(l - (1 - p)yv)r{y) (1 - yv){l - (1 - 2p)yp, - (1 - p)yu) 



%sinh(^ 



2y 2 ( e » — 1 



2y 3 I e v — 1 



(i/y)° 



(El) 

The previous expression for E 2 {p, v) has three terms, that we will call respectively A, B and C. We begin with A. 
Using the expansion 



sinh ; 



^ 2(l-2^- 1 )B 2j - __ 2 , 



(2i)! 



and the fact that 



we find 



r(y) 



.4 



(- 

y 



(\- p)v 2 (p+(l- p)v) 
2( e v+(i-p)v - l) 



(E2) 



(E3) 



(E4) 



We now calculate B. Using the expansion and recalling that all the odd Bj arc equal to except B\ = —1/2, we 
see that 



B 



(i-p)" 2 ( + p + (i-p)i 

2 2 



- (i - p)yv)r(y) 

4y 2 



(E5) 



We used here (|E3|) once again. We need the expansion of r(y) when y — > oo (again, after the expansion near // = 
and i/ = 0). At each order in p and v, r(y) is a polynomial in y. We have 



r(y) = 1 - ((1 - 2p)p + (1 - + 2p(l - p)(p 2 + pv)y 2 + O (y 3 ) 



(E6) 
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which is to be understood after the expansion in powers of p and v as usual. It finally gives 

2 (eM+Ci-rt" - 1) 2 ■ 1 U 

Thus, A and $ only contribute — p(l — p)p 2 v/2 to the eigenvalue. We will now see that C has a non-trivial contribution. 
Defining 

f(y) = —A — . (E8) 

e v — 1 

the term C of equation (|E1[) can be written 

r ^((1 - (1 - g>Q ,,, ^ (1 - 2p)p + v + (1 - p> 1 

C = g L/WV) 2 [/(y)](y2) + ^ L/(2/)](y3) • ( E9 ) 

Expanding /(y) in powers of r(y) with (|56[) . and expanding the powers of r(y) in powers of y, we obtain 

M = EEE |f ( k 7 *) 0) ("l)^-^^! - 2p)„ + (1 - (1 - p)vfi- 2 V j+1 - k - 1 ■ (E10) 

We take the term y r in the previous equation, setting i to 2j + 1 — k — r provided that it is nonncgative, and we move 
out of the sum over k the term k = 1 which is the only odd k such that ^ 0. We have 

[/<*, = J -f + f gw^ ('7)( 2j+1 ! a . r ) (Em 

x (-l) , '2« +1 - 2 *-'|(l -2p),i+(l -rt.f' +1 - 2, -'[,, + (l -p)^]"+ 2 ''- 2 >- 2 . 
For r > 0, j > fc + implies j > k (j integer). Thus 



k-i\ _ (-lp+ fe (2fc)!(2j-2fc)! 



j J 2 2 ij!fc!(j - fc)! 



(E12) 



Replacing j by j ; 4- k 7 we obtain 

oo 

E 



\ffM V- ^ (-lV +r (2j)! (fc-j)! / k \ 



2 2fc+r-i ( 2 j + l-r)! (k + r-j-l)l \j 
k=0 „-_T'--n w y v J y v 



x [(1 - 2p) M + (1 - p)^ +1 - r [p + (1 - p)v 



where [~(r — l)/2] is the smallest integer larger than (r — l)/2. For r = 1, the previous formula gives, rcsumming the 
sum over j 

lf(y)W) = E iS^ 1 - V + H fc ■ (EM) 



For r — 2, we obtain 



and for r = 3, we have 



klkV 

k=0 



lf(y)J(y*) = - E il^ 1 - V + ^) fc [(! - 2 ^ + C 1 - PM > (E15) 

fe=0 



oo _ OO „ 

[^)](v 3 ) = -E^^( 1 -^^ + Hl(l-2^+(l-p)^ + 2^T^^ • (E16) 

Inserting the last three equations into the expression (|E9[) for C, we see that everything cancels except the second 
term with r = 3. Thus, we have 

C = E^ff)iAl-^V+H fe . (E17) 
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Gathering everything, we finally obtain 

= + E fcfev^Cl - P) V + H* • (E18) 

fc=i '* 

This concludes the proof of the next-to-leading order of equation © . 
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